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ABSTRACT: We estimate the thermalization time in two colliding shock waves holo¬ 
graphic model of heavy-ion collisions. For this purpose we model the process by the 
Vaidya metric with a horizon defined by the trapped surface location. We consider two 
bottom-up AdS/QCD models that give, within the colliding shock waves approach, 
the dependence of multiplicity on the energy compatible with RHIC and LHC results. 
One model is a bottom-up AdS/QCD confining model and the other is related to an 
anisotropic thermalization. We estimate the thermalization time and show that increas¬ 
ing the confining potential decreases the thermalization time as well as an anisotropy 
accelerates the thermalization. 
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1 Introduction 


Holographic duality [1-3] provides a powerful tool for studying static properties of the 
QGP as well as its thermalization [4-6]. There are holographic models that reproduce 
perfectly the static properties of the QGP, meanwhile others holographic models are 
used to get non-static characteristics such as the thermalization time in heavy-ion 
collisions and the charged multiplicity. Holographic thermalization means a black hole 
formation in the dual space-time and particle multiplicities is defined by the entropy 
of the produced black hole. 

The gauge/string duality [1-3] perfectly works for the J\f—A Supersymmetric Yang- 
Mills theory, while a true dual description of real QCD is unknown, in spite of a lot 
of efforts have been made to find holographic QCD from string setup [7-13]. This ap¬ 
proach is known as the ”top-down” approach. Other approach, known as the ’’bottom- 
up” approach, is supposed to propose a holographic QCD model, i.e. a suitable vac¬ 
uum background, that fits experimental data and lattice results. The quark confining 
backgrounds that reproduce the Cornell potential, p-meson spectrum etc. have been 
proposed in [14-19]. Improved holographic QCD (IHQCD) that reproduce the QCD 
/3-function have been proposed in [20, 21]. The dual holographic approach has been 
successfully used to describe the static properties of the QGP, see [4] for review. 

The problem of the QGP formation in HIC is the subject of intensive study within 
the holographic approach in last years (see reviews [5, 6, 22] and refs therein) and 
initially has been considered in the AdS background [23]-[31] , that cannot describe 
either quark confinement either reproduce the QCD /3-function. However, as we have 
just mentioned, there are backgrounds which solve one, or even two of these problems. 
Therefore, to describe the holographic thermalization in more realistic frameworks, it 
is natural to study thermalization in these backgrounds. 

Thermalization in the improved holographic background has been studied in [32] 
and it has been shown that without additional assumptions, such as an energy depen¬ 
dent cut-off at the high energy [28], one cannot reproduce the multiplicity dependence 
on energy observed at RHIC and LHC in this background. In [33] it has been noticed 
that the holographic realization of the experimental dependence of multiplicity on the 
energy [35] requires an unstable background. In [36] it has been shown that the model 
reproducing the Cornell potential also gives an observed energy dependence of multi¬ 
plicities if one assumes that the multiplicity is related to the dual entropy produced 
during a limited time period. However in this consideration there is a limitation on the 
possible energy of colliding shock walls [36]. Since in this consideration we have used a 
more or less general background reproducing AdS at UV and confinement at IR, we can 
think that just general assumptions about the background prevent to reproduce suit- 
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able behavior in UV and IR and in the same time give the correct energy dependence 
of multiplicity at high energy. In particular, we can think that a default assumption 
of the isotropic form of background metrics is responsible for this discrepancy and for 
a more realistic description of entropy production in holographic models one has to 
consider anisotropic backgrounds [38]. 

In favor of anisotropic holographic backgrounds there are also additional arguments 
(see review [39] for holographic studies of strong coupled anisotropic theories). Up to 
a year ego, it was believed that just after heavy-ion collisions, a pre-equilibrium period 
exists for up to 1 frn/c and then the QGP appears and this QGP is isotropic. However 
now there is a belief that the QGP is created after a very short time after the collision, 
Ttherm ~ 0.1 fm/c, and it is anisotropic (’’anisotropic” means to a spatially anisotropy) 
for a short time r after the collision, 0 < Ttherm < r < Ti SO , and the time of locally 
isotropization is about Tj SO ~ 2 fm/c [37]. In the holographic version of this setup it 
is suitable to consider a black hole formation in a spatially anisotropic background. 
Motivated by recent experimental indications in favor of anisotropic thermalization, 
we also discuss a holographic thermalization scenario in the anisotropic 5-dimensional 
Lifshitz-likc background. The collision of domain walls in this background has been 
recently considered in [38], where it has been shown that for the critical exponent 
specifying the Lifshitz-likc background equal to 4, the dependence of multiplicity on 
the energy is desirable A 1 / 3 . 

In this paper we estimate thermalization time for two colliding shock waves in 
different backgrounds. Our main idea is very simple - the black hole creation in two 
shock waves collisions is modeled by Vaidya metric with a horizon corresponding to 
the location of the trapped surface appearing in two shock waves collision and thermal¬ 
ization time is estimated within standard prescription in the Vaidya metric [41]. The 
Vaidya deformation of isotropic backgrounds have been successfully used in description 
of thermalization in several isotropic models [41]-[50] as well as Vaidya deformations 
[51-54] of anisotropic metrics [56, 58-63]. 

As mentioned above, only for a special background the entropy of the black hole 
produced in the domain shock wave collision reproduces the energy dependence of 
particle multiplicities obtained at RHIC and LHC. We estimate the thermalization time 
for these cases. Namely, we estimate the isotropic/anisotropic thermalization time in a 
holographic bottom-up AdS/QCD confinement background that provides the Cornell 
potential and QCD /3-function. Here we use the Vaidya deformation of confinement 
background metrics in both isotropic and anisotropic cases. We also compare these 
results with our previous results obtained by general causality arguments [36] 

The paper is organized as follows. In Section 2 we discuss the isotropic case and 
in Section 3 the anisotropic one. Both sections are started by setup, where notations 
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and review of previous results are presented. In main parts of the sections estimations 
of thermalization time by the Vaidya modeling are performed. A comparison of results 
obtained by Vaidya modeling and by general causal arguments are presented in the end 
of these sections. 

2 Thermalization in isotropic backgrounds 

2.1 Setup 

2.1.1 General isotropic metric 

In a general isotropic holographic approach, the 5-dimensional metric is 

ds 2 = b 2 (z)(—dt 2 + dz 2 + dx 2 ), i — 1,2,3. (2.1) 

Following [32] we consider the following form of the b-factor 

P cz 2 /4 

*(*) = —, (2-2) 

where a and c are some constants. Metric (2.1) in the top-down approach supposes to 
solve the 5-dimensional dilaton-gravity equations of motion [20, 21]. It is also consid¬ 
ered in the bottom-up approach, in particular, the confining metric considered in [14] 
corresponds to a = 1 and c = 0.42 GeV 2 , see also [16-19]. 

2.1.2 Thermalization due to shock waves collision 

Point-like shock waves are usually considered in top-down backgrounds and they are 
supposed to be solutions of the 5-dimensional (dilaton) gravity with point-like sources. 
In the case of the point-like shock wave the deformation of metric (2.1) has the form 

dsthock = d s ~ + fids 2 shock , 

5ds 2 shock = b 2 (z)(t)(z ) x l ,x 2 )5{u)du 2 . (2.3) 

Here and below u,v — t±x 3 and the point shock profile cp(z, X \, x 2 ) solves the equation 

□ 3 60Os, Xi,x 2 ) = -16ttG 5 J uu , n 3b = (d a d a + 3^ + d 2 ), (2.4) 

here Cbb is the Beltrami-Laplace operator corresponding to metric ds\ b = b 2 (z)(dz 2 + 
dx 2 j_), _!_= 1, 2, G 5 is the 5-dimensional gravitational constant and J uu is the component 
of the energy-momentum tensor sourcing the point-like shock wave, J uu ~ 5(u)5(z — 
L)S(x_ l). 
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Two colliding point-like shock waves 

Sd 4-shocks = b2 ( z ) (0(*> x h x 2 )5(u)du 2 + (j)(z , x u x 2 )5(v)dv 2 ), u, v < 0, (2.5) 


produce a black hole, whose entropy S can be estimated by the area of the trapped 
surface 

5 > S TS , Sts = f A/ det I gsb\dzdx 2 ± = / b(z) 3 x]_(z)dz. (2.6) 

5 J z a h J z a 


z a and Zb are the points where ^{z a ,b) = 0 . ^(z) is the trapped surface profile function, 
which for the central collision [24, 32] up to the boundary conditions (that in fact define 
z a fi) satisfies the same equation as the shock wave profile. 

The simplest form of the black hole in the background metric (2.1) is given by [20] 


ds 2 = b~(z)(—f(z h ,z)dt 2 + 


dz 2 


f(z h ,z ) = 1 - K(z h ,z), K(z h ,z ) = 


f{z h ,z)z 2 
K{z) 


+ dx 2 ), i — 1, 2, 3, 


K( Zh y 


K(z) = 


dz 

b(zf' 


here Zh is the position of the horizon, and the temperature and entropy are 

1 47T 


T f\z h ) 
b 3 (z h )V 3 


= 47T 




b(z) ! 


S, = 


4 G, 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


where V 3 is the volume of the 3-dim space. 

For the shock domain walls [30, 31] the wave profile does not depend on the transver¬ 
sal coordinates and solves 


d l + Y dz 


(*) = 


16iyG 5 E S(z — z*) 


L 2 


b 3 (z) 


( 2 . 12 ) 


The trapped surface is located between points z a and Zb which satisfy equations [33] 


8nG 5 E 3 


(-a) 



8 vr G 5 E 
L 2 



b~ 3 dz 



(2.13) 
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here z* is the position of the collision point. The area density s (per the area in the 
transversal direction) of the trapped surface located between points z a and Zb, is given 
by 


1 


c Zb 


S = 


2 Gs 


b 3 dz. 


(2.14) 


Equations (2.13) give the relation between points z b , z a and the energy E, 

8ttG^E 


b\z a ) + b 6 (z b ) = 


L 2 


(2.15) 
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2.2 Thermalization Times 

2.2.1 Estimation with the Vaidya metric 

In this section we model the black hole creation in two shock waves collision by the 
Vaidya metric with a horizon corresponding to the location of the trapped surface 
appearing in these shock waves collision and estimate the thermalization time within 
the standard prescription in the Vaidya metric [41]. 

We relate z a and Zb, defining the location of the trapped surface, with the masses 
M a and Mb of the black brane in the background (2.1), 

M a = M(z a ), M{z a ) = K~\z a ), (2.16) 


where K(z) is defined by (2.9), and the same for Zb- We assume that from z a < Zb 
it follows that M b < M a . The appearance of the trapped surface located at z a and 
Zb = oo can be modeled by the Vaidya metric 

ds 2 = b 2 (z) (—/( Zh , z, v ) dv 2 — 2 dvdz + dx 2 ) , (2-17) 

f(z a ,z,v) = 1 -9(v)K(z a ,z). (2.18) 


As it is accepted in the Vaidya approach [41]-[50], to find the thermalization time 
r at the scale i one has to consider a geodesic with equal time endpoints located at the 
boundary at distance £ and find the time r, when this geodesic is covered by the black 
shell (2.17). The thermalization time r and the distance i are related as 


£ = 


T = 



1 b(s) dw 

b ^ SW ^ 1 - K(z h ,sw)) 


dw 

1 - K(z h , sw) ’ 



&2 (s) A 

b 2 (sw) J 


(2.19) 


( 2 . 20 ) 


i.e. one finds the thermalization time at the scale l excluding the axillary parameter s 
from the system of equations (2.19) and (2.20), 


r = 


Ttherm (^hi O' 


( 2 . 21 ) 


Assuming z a < Zb we estimate the thermalization time due to the trapped surface 
formation by 


7~therm(Zai 0 T~therm\Zai O' 


( 2 . 22 ) 
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2.2.2 Estimation for confining metric 

In this section we consider thermalization for the metric with the confining 6-factor 


b{z) = (2.23) 

here c is related with the notations of [15] as c = \caz and we use the bottom-up 
version of the blackening function considered in [14] 


K(z) = z 4 , f(z h ,z) = l--r. 


(2.24) 


Note that according (2.16) there are non-leading corrections to (2.24). Indeed, for the 
6 -factor (2.23) the blackening function according (2.16) is 


K(z a , z) = 


dz 

b(zY 


-l+e~ 3cz2 +3e~ 3cz2 cz 2 
-1 + e~ 3cZh2 + 3 e~ 3cZh - 2 czh 2 


(2.25) 


The leading term is in agreement with (2.24) 


K(z a , z) = k±z 4 + k(>z 6 + O (z 8 ) , (2.26) 

where k± = l/z£ + 2 c/z\ + O (1), k e = —2c/z^ + 0(1). 

In Fig.l.A the dependence of the thermalization time r on the scale i for the 
metric (2.23) for different values of c is shown. From this plot we see that increasing 
the confinement potential, in fact when c in formula (2.23) increases, we decrease the 
thermalization time. 

In Fig.l.B. the dependence of v — t/r on l for the same parameters as in Fig.l.A 
is presented. We can interpret v = l/t as a propagation velocity of the thermalization. 
We see that the velocity increases with increasing of the factor c in (2.23). 

In Fig. 2. A. dependencies of r on t for different masses of the shell and the same 
parameters c as in Fig.l.A are presented. We see that for chosen parameters the 
dependence on Zh is very small. We see that increasing Zh for AdS case we increase the 
thermalization time. The same is true for c = 0.1. However, when we increase c so that 
0.2 < c < 0.5 the dependence on Zh becomes more essential (the distances between the 
magenta lines are larger as compared with distances of green lines) and when increasing 
Zh. we decrease the thermalization time. We also see, Fig.2.B., that in considered cases 
only up to some given distance the thermalization is possible. This is related with the 
breaking of geodesics with two large endpoints distance in the confining background. 





rL fm] v 



Figure 1 . A. Dependencies of t on l for 5-dimensional metric (2.17) with the confining 
6 -factor (2.23) for c = 0 (red), c = 0.1 (blue), c = 0.2 (green), c = 0.5 (magenta), c = 2.56 
(cyan) , c = 5.16 (brown ) and the blackening factor (2.16) with Zh = 1. B. Dependencies of 
v = 1/t on i for the same parameters as in the left panel. 
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Figure 2. A. Dependencies of r on l for the 5-dimensional metric (2.17) with the confining 
6 -factor (2.23) for c = 0 (red), c = 0.1/m ~ 2 (blue), c = 0.2/m -2 (green), c = 0.5 fm~ 2 
(magenta), c = 2.56 /m -2 (cyan), c = 5.16 /m ~ 2 (brown ) and Zh = 1/m (solid lines), 
Zh = 1.2/m (dotted lines), Zh = 1.8 fm (dashed lines). B. The zoom of the plot A in the 
region 2 fm < £ < 9 fm and 1 fm < r < 1.8 fm. C. The zoom of the plot A in the region 
0 < £ < 8 fm and 0.3 fm < r < 0.5 fm. 


In Fig. 3 the thermalization process is presented for more realistic values of c, namely 
c = 2.56 f m~ 2 . Dependencies of r on t for 5-dimensional metric (2.17) with the confin¬ 
ing 6 -factor (2.23) with c = 2.56 and different blackening functions: blackening function 
(2.28) for Zh = 1.2 and Zh = 1.8 corresponds to solid and dashed blues lines; blackening 
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Figure 3. A. Dependencies of r on £ for the 5-dimensional metric (2.17) with the confining 
6-factor (2.23) with c = 2.56 and different blackening functions: blackening function (2.28) 
for Zh = 1.2 and Zh = 1-8 corresponds to solid and dashed blue lines; blackening function 
(2.25) for Zh = 1.2 and Zh = 1.8 corresponds to solid and dashed red lines. B. The zoom of 
the plot A. C. Dependencies of r j£ on l for the same parameters as in the plot A. 


function (2.25) for Zh = 1.2 and Zh = 1.8 corresponds to solid and dashed red lines. 
We see that corrections (2.26) do not play an essential role for small distances, but at 
large distances the blackening factor (2.24) admits longer geodesics as compared with 
(2.25). 


2.2.3 Estimation for power metric 

It is instructive to compare estimations obtained in the previous subsection 2.2.2 with 
results obtained for the intermediate metric [36] 



In this case the blackening function is 

~da -\-1 

/£> 

„dci +1 ’ 

Z h 



(2.27) 


(2.28) 


and the thermalization time r 


£ = 2 s [ 

Jo 



the scale l is 

w a dw 

dw 

(1 - ff£w 1+da ) 

ft 


(2.29) 

(2.30) 
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Figure 4. A. Dependencies of r on i for 5-dimensional metric (2.17) with the power-law 
6 -factor (2.27) for a = 1 (solid red line), i.e. the AdS^ case and a = 0.5 (thick gray line). B. 
Velocities for the same parameters. 


In Fig.4.A dependencies of r on l for 5-dimensional metric (2.17) with the power 
6 -factors (2.27): a = 1 (solid line), i.e. the AdS 5 case, and a = 0.5 (dashed line) are 
presented. We see that decreasing a we decrease the thermalization time. From Fig. 4 
we see that the velocity of propagation of thermalization increases with decreasing a. 

It is interesting to note that the dependence of t on l does not depend on the 
position of the horizon for a — 1 , i.e. 


’■(zk.Olo-l = t(1, Ola-1 


(2.31) 


The confirmation of (2.31) is presented in Fig.5. Indeed, we see in the plots in 
Fig.5 that three lines, the purple thick line, the red solid line and the dashed blue one 
coincide, i.e. the thermalization time does not depend on Zh- 
From (2.29) and (2.30) it is evident that 


— = 2S 

Zh 


^ = S 


w a dw 


-\/(l — S 1+da w 1+da ) (1 — w 2a ) 
dw 


(2.32) 


(2.33) 


Zh J Q (1 — gi+da w l+da^ ’ 

here we introduce a new parameter S = s/zh- The system of equations (2.32) and 
(2.33) is nothing but the system of equations that defines the thermalization time in 
units of Zh for the Vaidya model with the unique mass, i.e. 


T~ Zh Z"therm (1 1 _ )• 


(2.34) 


Note that in [36] we have shown that the b-factor for the intermediate metric 
(2.27) with a = 0.5 and L e ff = 20.86 approximates the confining b-factor for 1.2 fm < 
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Figure 5. A. Dependencies of r on i for 5-dimensional metric (2.17) with the power 6-factor 
(2.27) for a = 0.5 and blackening factor (2.28) with different Zh : Zh = 1 (blue dashed line), 
Zh = 1.5 (red solid line), Zh = 2.5 (purple thick solid line). 



Figure 6. A. Dependencies of r on i for 5-dimensional metric (2.17) with different 6-factors: 
the thin red line corresponds to power factor (2.27) with a = 1 (AdS); the gray line corresponds 
to a = 0.5 and the thick cyan line corresponds to the factor (2.23) with c = 2.56. B. 
Dependencies of velocities of the propagation of thermalization on the characteristic distance. 

z < 1.8/m. For z ~ 1.1 fm the thermalization times for these two models are very 
close, however, for the distance z ~ 2 fm thermalization times are different up to 
factor 0.3. For larger distances these two models show essentially different behaviour: 
the thermalization time for the confining metric does not depend essentially on the 
distances, meanwhile the intermediate model has approximately a linear dependence 
of the thermalization time on the distance. 
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In Fig. 6 . A. we plot the dependence of the thermalization time for 3 different 
models: the AdS case (thin solid red line), the intermediate case [36] (gray line) and 
for the confining factor [14] (thick cyan line). For all cases z b — 1 /m. 

The small dependence of the thermalization time on the value of the mass of the 
shell for special models can be seen from the following considerations. Assuming that 
b(z ) satisfies the scaling 


b[kz) = k^b^z) + ..., 0 < z < Zh, (2.35) 

one can reduce the problem of finding the thermalization time at the scale l by the 
Vaidya metric with mass Mh to the same problem with the unite mass. Indeed, from 
(2.35) and (2.16) the scaling for the mass M(z) is 

« k^-'M-Hz). (2.36) 

k 

Performing the rescaling 

z = k a z, x = k a x, t = k a t, (2.37) 

with 

k a = Mt^ X 

we recast the metric (2.17) to the same form in terms of the 
blackening function with a unit mass f(l,z,v), 



b / „ \ 
ds 2 = Mt b 1 b 2 (z) (—/(1, 5, v) dv 2 — 2 dvdz + dx 2 ) . 

(2.39) 

Therefore, 

'^'therrrS^'ai^) ^ aether ^aTthermi^-') ^/^a) • 

(2.40) 

Assuming that 

Ttherm{ lj^) = t + C' 2 ^ + ••• 

(2.41) 

we get 

T t herm(za, Z b , £) ~ Ci + C 2 f/k 2 a . 

(2.42) 


2.2.4 Estimation of the formation time of a trapped surface by non-local 
correlation functions 

We can also estimate the trapped surface formation time by a characteristic size of the 
trapped surface [36], i.e. 

T,her m ~ ~— • (2.43) 

V z 


(2.38) 

tilde-coordinates and the 
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Here v z is the velocity of propagation of a signal along the ^-direction. We can estimate 
the velocity of propagation of the signal along the z direction in two ways. We can 
relate two points on the boundary by the geodesic or the string worldsheet stretched 
on the static quarks world lines located at these points. The geodesic as well as the 
worldsheet have maximum z-coordinates, z*. Varying the end points we vary z* and 
we can estimate 


v = 


A z* 

Ax 


■ c, 


(2.44) 


here we take into account that any propagation in the ^-direction is limited by the light 
velocity. Below we assume c — 1. 

Estimation with string. The relation between the interquark distance x and the 
string maximum holographic coordinate z m is given by 


From this formula we have 



(2.45) 


A^ _ / b\z) 

Ax V b 4 (z m ) 


(2.46) 


We see that (2.46) defines the z-dependent velocity. In particular, considering this 
estimation in the confining background [36, 40] where 1.2 fm < z m < 1.8 fm and 
we get Ax = Az m /2A. This estimation gives the trapped surface formation time 
Ttherm ~ 0.25 flTl. 

The estimation with geodesics in the intermediate background gives Ax = |A z m , 
which is 4 times longer as compared with the string estimation. 


3 Thermalization in anisotropic backgrounds 

3.1 Setup 

3.1.1 Anisotropic metrics 

As an anisotropic background (with the space anisotropy) we consider a five-dimensional 
Lifshitz-like metric [56, 57, 62] 


ds 2 Ll = L 2 


(— dt 2 + dx 2 ) ( dy\ + dy 2 ) dz 21 

z 2 ^ z 2 ! v + z 2 


(3.1) 


Let us remind that the anisotropic Lifshitz metric with space symmetry is given by 
[55, 56, 58-63] 

-dt 2 (dx 2 + dy\ + dy 2 ) dz 2 


ds\if — L~ 


+ 


^ 2 /^ 


+ 


(3.2) 


14 
















Both metrics for v = 1 are reduced to the Poincare patch of AdS 5 . 


3.1.2 Thermalization due to shock waves collision 

This scenario assumes that the main part of multiplicity is produced in an anisotropic 
regime and this part of multiplicity can be estimated by the trapped surface produced 
under a collision of the two shock waves in an anisotropic background. This scenario 
is accepted in the recent paper [38], where collisions of shock waves in the Lifshitz-like 
background have been considered. 

As a model of anisotropic background we consider a five-dimensional Lifshitz-like 
metric (3.1). The shock domain wall moving in the u-direction is given by deformation 
of the metric (3.1) 

1 2 _ 7 2 , T 2 <KZ)S(U), 2 (o q\ 

d s DW,Li ~ d s Li + L du , (3.3) 

with the profile function cj)(z) satisfying 


d 2 4>(z) 
dz 2 



1 d(p(z) 
z dz 


/ z \ 1 + 2 / zy 

\GttG^J uu ^ J uu E ( — J d(z %*)•> 


(3.4) 


where z* is the ^-coordinate of the collision point. The entropy can be written in terms 
of z a and Zb defining the location of the trapped surface 


v 


s = 


4 G 5 \(z a y/" (,z b ) 2 A 


The analog of relation (2.15) is 


„-l-2/v , - 1 - 2 /" _ 

*a ' “b 


87 tGzE 


12/v+3 ■ 

The maximal entropy is achieved at Zb —>■ 00 and it is 

” -(8nG 5 ) 2/{u+2) E 2 ^ u+2 \ 


s = 


4Gb 


(3.5) 


(3.6) 


(3.7) 


The leading asymptotic gives rise to the value of multiplicity, which is the most 
compatible to the experimental data, for v = 4 [38]. Note, that as in the isotropic case 
for intermediate values of the energy we have to take into account the next to leading 
term, but now we do not have a restriction on z a from below since we do not obliged 
to fit our metric to the metric with a given b- factor. We have changed the background 
on which we consider the collisions of domain walls. It would be interesting to find 
’’top-down” motivations for consideration the anisotropic background as a holographic 
model for HIC. 
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3.2 Thermalization time in the anisotropic background 
3.2.1 Estimation with Vaidya metric 

We can estimate the thermalization time in the anisotropic background 


ds 2 
f(z h ,z ) 


; 2/ \ ( f ( Z hi z) , 2 

6 (*) [~^=w dv 
1 g(g) 

K(Za) ' 


dvdz 


yV— 1 



For metric (3.8) with f(zh, z) as in (3.9) we have 


(3.8) 

(3.9) 


1 = 



b(s) 

b(sw) 


dw 



T = 


' 0 W 


1-v 


dw 



K(sw) \ 


(3.10) 


(3.11) 


For the power-law 6-factor, b(z) = ( L/z) a the blackening factor is given by (3.9) 
with 

K(z) = z ad+u (3.12) 

Thermalization for these models have been considered in [51, 52], In Fig.7 we show 
the dependence of the thermalization time on the distance for the model with a = 0.5 
and different Zh- We see that increasing the anisotropy the dependence on the horizon 
position increases as well. We also see that the breaking of the geodesics appears at 
smaller distances for the case of small Zh (heavy black masses). 

For the confining 6-factor by the analogy with (2.24) we use the ’’phenomenological” 
blackening function 

K(z) = z 3+1/ . (3.13) 

In Fig. 9 the influence of the anisotropy on the thermalization time for the confining 
metric (3.8) for different c is presented. We see that anisotropy essentially decreases the 
thermalization time. Meanwhile there is no essential dependence of the thermalization 
time on the mass of the shell. 


3.2.2 Estimation with non-local correlation functions 

To estimate the trapped surface formation time we use an analog of estimations (2.43) 
and (2.44). Considering the relation between the interquark distance along yj-direction 
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r[/m] 



Figure 7. A. Dependencies of t on l for 5-dimensional case for a = 0.5 and different v and 
different Zh : solid lines Zh = 1 and dotted lines Zh = 2. Blue lines: v = 2; magenta lines: 
zz = 3; purple lines: zz = 4. B. Dependencies of velocity of thermalization propagation on the 
distance. The zoom to the plot B. 


r[fin\ 



Figure 8. Dependencies of r on i for the 5-dimensional metric (3.8) and the blackening 
factor (3.13) for different values of c and v. Red lines correspond to the AdS case, c = 0 and 
cyan lines to the confining b-factor with c = 2.56 fm~ 2 . The isotropic case, v = 1, is shown 
by solid lines, anisotropic cases with v = 2 are shown by dashed lines, v = 3 by dotted lines 
and v = 4 by dot-dashed lines. For all cases z a = 1.2/m. 
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Figure 9. Dependencies of r on £ for 5-dimentional metric (3.8) and blackening factor 
(3.13) for different c, different is and different z a . Red lines correspond to the AdS case, 
c = 0: isotropic case (is = 1) solid z a = 1.2 and z a = 1.8 dashed lines, anisotropic case (is = 3) 
z a = 1.2 dotted and z a = 1.8 dot-dashed lines. Cyan lines correspond to the confining 6-factor 
with c = 2.65: isotropic case (is = 1) z a = 1.2 solid and z a = 1.8 dashed lines, anisotropic 
case (is = 3) z a = 1.2 dotted and z a = 1.8 dot-dashed lines. Brown lines correspond to the 
confining 6-factor with c = 5.75 and different anisotropy factors and different z a '- isotropic case 
z a = 1.2 solid and z a = 1.8 dashed lines, anisotropic case (is = 3) z a = 1.2 dotted and z a = 1.8 
dot-dashed lines. B. The zoom of plot A, 5.8 fm < £ < 6.8 fm, 0.027 fm < r < 0.028 fm. 


and the maximum of the string profile holographic coordinate z m we get the relation 


Vi = 



dy 

\/l — z 2 ~ 2 / v 


(3.14) 


In particular for v = 4 we have y\ = 2.22 z^ and therefore, Ay & 0.55 A ^J2\ that 
for the trapped surface located at the intermediate zone 1.3 fm < z < 1.8 fm gives 
'Tt.hermn ~ 0.2 flTlf C. 

Note, that the estimation with geodesics, gives y\ = 2 iszm, u , which is once again in 
four times longer as compared with the string estimation. 

One can compare these estimations of the trapped surface formation time with the 
thermalization time at scale about 2 fm -j- 4 /m, that is according sect. 3.2.1 for v = 4 
is ~ 0.05 fm. 
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4 Conclusion 


As it has been mentioned in the Introduction, the entropy of the black hole produced 
in the shock wave collision predicts multiplicities for heavy ion collisions at RHIC and 
LHC only for special backgrounds. In this paper we have estimated the thermalization 
time for these cases. 

In particular, we have estimated the anisotropic thermalization time in a holo¬ 
graphic bottom-up AdS/QCD confinement backgrounds that provides the Cornell po¬ 
tential and QCD /3-function. We have shown that thermalization time is up to 5 times 
faster comparing to the isotropic case. It is interesting that we have not seen essential 
dependence of the thermalization time on the temperature, i.e. our method predicts 
the same order of anisotropic thermalization time for RHIC and LHC. 
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